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Abstract 

The Stokes and anti-Stokes components in the inelastic scattering of light are related to phonon 
statistics and have been broadly used to measure temperature and phonon lifetimes in different 
materials. However, correlation between the components are expected to change the Stokes/anti- 
Stokes intensity ratio, imposing corrections to the broadly used Bose-Einstein statistics. Here the 
excitation power dependence of these scattering processes is theoretically described by an effective 
Hamiltonian that includes correlation between the Stokes and the anti-Stokes events. The model 
is used to fit available experimental results in three-dimensional diamond and two-dimensional 
graphene, showing that the phenomenon can significantly increase in the low-dimensional system 
under specific resonance conditions. By setting the scientific basis for the Stokes-anti-Stokes cor¬ 
related phenomenon, the use of the Bose-Einstein population function for reasoning the inelastic 
scattering is generalized, providing a model to predict the conversion of optical phonons into heat 
or light, according to coupling constants and decay rates. The model applies to inelastic scattering 
in general. 
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The inelastic scattering of light by matter [I] exhibits two components: the Stokes (S), 
where an incident photon is converted into a phonon and a red-shifted S photon, and the anti- 
Stokes ( aS ), where one incident photon and one existing phonon are annihilated, generating a 
blue-shifted aS photon. The Stokes/anti-Stokes intensity ratio is a signature of the quantum 
character of lattice vibrations, considered to be defined by the Bose-Einstein distribution 

l aS _ n n 0 

Is 1 + no 

where C depends on the optical setup, and no = (e fe s T — l ) -1 is the effective phonon 
population [2j. Hq 1 for phonons with energies (hu) larger than the thermal energy ksT 
{ks is the Boltzmann constant). The I a s/Is intensity ratio is used to measure phonon 
lifetimes [ 3 ], local effective temperatures IH El and optical resonances [6], and except for 
very specific resonance conditions, the I a s/Is ratio should approach the unit only if the 
temperature is high enough to activate a very large phonon population. However, Klyshko 
E !H| has proposed a correlated process, called here the Stokes-anti-Stokes ( SaS) event, 
where a phonon created by the Stokes process is subsequently annihilated in the anti-Stokes 
process. If the SaS event is significant, the picture described by the Bose-Einstein phonon 
distribution is not complete, and Eq.[T]has to be generalized. 

Evidences for the SaS process are accumulating in materials science MEHm, generating 
interests in quantum optics j9j lH]. Lee et al. [9] demonstrated the quantum nature of 
the SaS correlation in diamond by measuring a non-classical SaS field correlation function 
g^ 2 \ Klyshko [8] pointed that the correlated character of this S and aS photons can be 
continuously varied from purely quantum to purely classical and Kasperczyk et al. HU 
explored this transition in character by changing the photon and phonon reservoirs through 
changes in the excitation laser power. While these experiments are usually performed with 
ultra-fast pulsed lasers to enhance the response of the non-linear SaS event, Jorio et al. )10j 
have provided evidences for the observation of dominant SaS event using a few miliwatts 
continuum wave (CW) laser, i.e. achievable even with a simple laser pointer. This result was 
obtained in twisted-bilayer graphene (tBLG), a two-dimensional system specially engineered 
to exhibit resonance with the aS photon emission HDj. Since phonons have a significant 
lifetime, it has been proposed that these systems can work as a solid state quantum memory, 
storing information between the write ( Stokes ) and read ( anti-Stokes ) processes [9|. In 
diamond and graphene Tw fcgT, and the quantum memory would be able to work at 
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room temperature. 

Here the I a s/Is~ based phonon population analysis is generalized by proposing an effective 
Hamiltonian that explicitly considers the Stokes and anti-Stokes fields correlation, on a 
simple formalism that is able to fit the experimental results from diamond and graphene. 
The Hamiltonian is given by 

H = huodJa + hisAc + ftidsogbs + hid a sb' aS b a s 

+h\s(ac^b^ s + h.c .) + h\ aS [acb\ s + h.c. ), ( 2 ) 

where b s , aS (b aS ), c (cf) and a (cd) represent the annihilation (creation) operator of Stokes 
(S), anti-Stokes (aS), phonon (c) and incident ( a ) fields respectively. A5 and A a s are the 
coupling constants. The A5 coupling term describes the creation of an S photon and a phonon 
through the absorption of an incident photon, and the A a s term describes the creation of 
an aS photon through the absorption of an incident photon and a phonon, ojq and is are 
the incident and the phonon fields frequencies, respectively, then the Stokes and anti-Stokes 
modes have energies given by hu a s,s — bw 0 ± his. Typically, the coherent sources used in 
experiments of this type can be considered to have large enough number of photons, which 
allows us to replace (a, at) —>■ |a|, with |a| J being the mean number of incident photons. 
The laser power is given by Pl = Al|a| 2 , A being a constant with power units, depending 
on the laser frequency. The Hamiltonian model above is valid within the coherence time of 
the pumping laser, whether continuum or pulsed. 

Hamiltonian (| 2 | accounts for the interaction at the material but it does not properly 
describe the dissipation of the created excitations due to the presence of the photonic and 
phononic reservoirs. Such dynamics can be computed by means of the Markovian master 
equation H2| for the overall reduced density operator of the three fields, introducing decay 
rates of phonons and scattered photons. The master equation for the density operator 
p = Tr fl p to tai, in the Lindblad form, reads 

j t P = +C(p) ( 3 ) 

where Tr n indicates the tracing out of the reservoir degrees of freedom. The Lindbladian 
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term in our model is C = Cb + £ c , with: 


A,(p) = - 7s(6£M + p6£6 x «- 26 x/ 56t) 

x=S,aS 

£c(p) = ~lc{ n 0 + l)(c t cp + pc f c - 2cpc f ) 

— 7c'^o(cc^p + pcc* — 2c^pc), (4) 


with 75 , 7 a 5 and 7 C being the decay rates (proportional to the inverse of the coherence 
time) of the respective Stokes, anti-Stokes and phonon fields. To analyze the Stokes and 
anti-Stokes held intensities and their respective correlation function at zero delay (r = 0), 
we compute the average values 


(ns,os) = (b\,o,sis,o,s), (n c ) = (c'c), g 2 ( 0 ) 


(Kblsksh) 

(bibs) (Ksks)' 


( 5 ) 


The power dependence for the intensity ratio I a s/Is and the SaS held correlation g 2 ( 0 ) 
are described in panels (a) and (b) of Fig.[TJ respectively, for different values of the thermal 
phonon population no- While the Stokes held is found to be proportional to the excitation 
laser power Pl , the anti-Stokes held exhibits two different regimes. For lower power, ther¬ 
mal phonons dominate the process and (■ n a s) is proportional to Pl (assuming there is no 
laser induced heating). For higher powers, the SaS phenomenon dominates, and (n a s) is 
proportional to P 2 (see Supplementary Information for details). This rationale explains the 
results in Fig.[l](a), I a s/Is being constant for low power, and for high power, hs/Is oc Pl- 
In Fig.JlJb) , the SaS correlation function g 2 ( 0) is shown to be proportional to the inverse 
laser power, and it goes to a no-dependent constant value for small values of Pl- 

The behaviors of Is , I a s an d g 2 ( 0) observed by Kasperczyk et al. [llj are, therefore, 
perfectly described by the steady-state number of S, aS photons (■ ns ) ss , {n a s)ss an d g 2 ( 0 ) 
in Eqs.[5j calculated by considering the model in Eqs.[2]{5] Here, (■ ns >a s)ss is the steady state 
number of S,aS photons. 

A rationale for I a s/Is that takes into account the SaS phenomenon can then be developed 
based on the theory for ( n aS ) ss /(n s )ss■ Setting A s = X aS = X , 7 s = 7 as = 7, M a \hc ^ 0, 
and considering the limit of 7 / 7 c no, the anti-Stokes/Stokes intensity ratio is given by 
(see Supplemental Information) 


("as) "a 1 + Alt 

(,1 S ) n „ + l 1-55^1,^ 


( 6 ) 
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where Pq ~ M 77 C / 2 A 2 for 7/70 no, given in power units. In the experimental range of 
interest, Eq.[ 6 ]can be approximated by 



where the fundamental constants X 2 /A / y'y c = Cs a s measure the importance of the SaS 
phenomena, and C' = C a s/Cs, with Cg.aS dehned as the proportionality constant such that 
Is,a,s = Cs,as(ns,as)ss■ The constant C depends on the optical parameters, for instance, 
absorption coefficients, inelastic scattering cross sections and the dimension (geometry) of 
the sample under study m 

Figure|2] shows a comparison of the experimental I a s/Is results from the three-dimensional 
diamond (black data) | 11 | . from a two-dimensional bi-layer graphene, where the layers are 
superposed in the so-called AB-staking configuration (AB-BLG, open red data) [12], and 
finally from the two-dimensional twisted bi-layer graphene (tBLG, solid red data), which 
was engineered to exhibit electronic resonance with the aS photon emission [10 ], 

The three different results observed for the three different samples in Fig. [2] can all be 
fit with Eq.[7j giving the parameters provided in table]!] (see Supplementary Information 
for fitting details). Analysis of the data shows that the increase in I a s/Is f° r diamond is 
dominated by the SaS process ( Cs a s — 2.13 x 10 -5 mW -1 ). Actually, the diamond data 
can be satisfactorily fit by imposing a constant temperature ( T = 295 K), with a larger SaS 
contribution ( CsaS — 3.9 x 10 -5 mW -1 ), consistent with the fact that, for measurements 
with a CW laser, within the same excitation power range, the I a s/Is power dependence in 
diamond exhibits a constant behavior (see data in the Supplementary Information). For 
the AB-stacked bilayer graphene (AB-BLG), the I a s/Is power dependence is dominated 
by heating, with a linear dependence of the effective phonon temperature on laser power 
(T[K] = 295 + 30Pi[mW]). Finally, the increase in I a s/Is f° r twisted-bilayer graphene 
(tBLG) is dominated by the SaS process ( Cs a s — 1-26 x 10 -3 mW -1 ), with a contribution 
from laser induced heating (T[K] = 295 + 19 Pl[iuW]) that is consistently smaller than that 
from AB-BLG. Notice that the constant found to rule the temperature increase in tBLG 
(19 ± 2K/mW) is roughly half of the value obtained for the AB-BLG (30.1 ± 0.2K/mW), 
consistent with the G band frequency changes measured in tBLG and AB-BLG [TO]. Due 
to resonance with the aS photon emission in tBLG, the assumption A 5 = A a s should be 
discussed. However, considering A a s 7 ^ A 5 does not change the overall picture. More details 


5 







can be found in the Supplemental Information, including a comparative analysis of the 
thermal versus SaS contributions to I a s/Is f° r the three samples. 

Comparing the graphene samples, while in AB-BLG the generated phonons are mostly 
converted into heat, in tBLG they are mostly converted into light (aS photons). But the 
most striking result is the much larger efficiency of the SaS process in tBLG, with a CsaS — 
(A 2 / 77 c )A~ l value that is two orders of magnitude larger than in bulk diamond. Besides 
these two orders of magnitude difference, the Pf dependence for I a s in diamond can only 
be observed with the use of femtosecond lasers, where the pulse intensity is about 10 5 times 
larger than the intensity of continuous wave laser radiation of the same average power. 
Furthermore, in bulk diamond the number of active atoms in the focal volume is ~ 10' 
times larger than in the two-dimensional graphene system. Finally, the ratio I a s/Is at the 
highest power in tBLG reaches values of ~0.5, while in diamond I a s/Is ~ 0-1 was reached. 
Putting all numbers together, correlated SaS generation per involved atom in the graphene- 
based system is roughly 10 14 times more efficient than in bulk diamond. 

Qualitatively, the reason for the striking efficiency of correlated SaS Raman scattering 
in the twisted bilayer graphene system might be related to a reduction in phase-space in 
the scattering event. In real space, the confined two-dimensional structure, which enhances 
electron-hole interactions in low-dimensional structures dans], leads to a large overlap be¬ 
tween the photon and phonon wavefunctions. In reciprocal space, the logarithmically diverg¬ 
ing two-dimensional Van-Hove singularity at the hexagonal saddle (M) point |16| strongly 
enhances the generation and recombination of electron-hole pairs at the aS photon energy 
[nmm. The SaS event is also expected to play a role in other low dimensional materials 
[18], where phase-space reduction becomes important. 

Raman spectroscopy is established as an important tool to study and characterize nanos¬ 
tructures [2j. By setting the scientific basis for the SaS phenomenon, the use of the Bose- 
Einstein population function for reasoning the inelastic scattering is generalized, providing a 
model to predict the conversion of optical phonons into heat or light, according to coupling 
constants and decay rates. It is surprising that the correlation between the Stokes and the 
anti-Stokes components can seriously affect I a s/Is, being measurable indirectly from the in¬ 
tensity ratio. This aspect made it possible to measure the importance of the SaS correlation 
in graphene, where a direct g 2 ( 0) measurement is not possible due to ultrafast luminescence 

mm- 
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The stationary phonon population can then be analytically obtained from the steady 


state solution of Eq.(|3| and, in the limit of 7/y c n 0 , it is given by 



77 c 


( 8 ) 


The phonon population has two contributions, one thermal and the other related to the 
SaS phenomenon, each contribution dominating at different excitation power ranges (see 
a plot of (n c ) ss in Fig. Is of the Supplementary Information). The coupling constants and 
the decay rates governing the phenomena can themselves be dependent on temperature, 
through anharmonicities or electron-phonon scattering. On the other way around, the SaS 
event is an important decay channel that has to be explicitly considered when using Raman 
spectroscopy to extract structural and transport properties related to phonon scattering in 
low-dimensional structures [2Tb23l| . 

More generally, since the Stokes and anti-Stokes components are common to inelastic 
scattering processes in general, the model is not limited to the Raman scattering by phonons. 
It is generic, for instance applying also to the scattering of He Atoms on surfaces, to the 
inelastic neutron scattering, etc. Experimental and theoretical work is needed to address 
the physics of the A s, a s, 7s ,as and y c parameters when applied to different materials and 
scattering phenomena. 
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Table I: Parameters used to fit the data from the three materials depicted in Fig. [2] using Eq.[7] 
CsaS is given in mW _1 (C', A, 7 and r y c are dimensionless, while A is given in mW). Temperature 
T is given in K and laser power ( Pi ) in mW. The constant multiplying Pi is given in K/mW. 


MATERIAL 

a 

CsaS = (A 2 / 77 c )M -1 [mW -1 ] 

T[K] 

diamond 

(10 ± 1 ) 

(2.13 ± 1.04) x 10 " 5 

295 ± (0.37 ±0.15)P L 

tBLG 

(20 ± 1 ) 

(1.26 ± 0 . 01 ) x 10" 3 

295 ± (19 ±2 )P L 

AB-BLG 

(3.5 ±0.2) 

(0.24 ± 1.36) x 10 -5 

295 ±(30.1 ±0.2 )P L 
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Figure 1: Theoretical description of the Stokes and anti-Stokes intensity and correlation phenomena, 
(a) gives the population ratio (n a s) S s/( n s)ssi where ss stands for steady state, and (b) gives the 
SaS field correlation g 2 { 0), according with Eqs.[5j The excitation laser power Pl dependences 
are plot for different values of the thermal phonon population no (see legend in (b)). We set 
A s = A as = A, 7 5 = 7 aS = 7 and Pq « 7 l 77 c / 2 A 2 for 7/70 no- Po in power units. The stationary 
phonon population gives nonphysical description for Pl > Po (more details in the Supplemental 
Information). 
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Figure 2: Fitting the I a s/Is intensity ratio in three different solid state systems: three-dimensional 
diamond (black data), two-dimensional bi-layer graphene (AB-BLG, open red data), and two- 
dimensional twisted bi-layer graphene (tBLG, solid red data). Data points are experimental results 
from Refs. [TT] (diamond) and |10| (graphene-systems), and lines are fit to the data using Eq.[7]and 
the parameters listed in table]]] 

SUPPLEMENTAL INFORMATION: POWER DEPENDENCE OF KLYSHKO’S 
STOKES-ANTI-STOKES CORRELATION IN THE INELASTIC SCATTERING OF 
LIGHT 

Here we provide more details about the theoretical model and include some results not 
shown in the main text. We also present a discussion about the analysis of the experimental 
data from the literature. 


MORE DETAILS ABOUT THE MODEL 


To analyze the Stokes and anti-Stokes populations and their respective correlation func¬ 
tions at zero time delay (r = 0), we computed the average values 


(■ n s(aS )> = (bl {aS) b S (aS)), (™c) = (^c), g\ 0) = ’ 

( b S b s)( b aS b aS) 


(9) 
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which can be obtained by means of the equation 


j t (6) = i([H,6]) + M6m, ( 10 ) 

yielding the following linear system of differential equations: 

d t {n s ) = -2X a lm{(cb s )} - 2^ s (n s ) 
d t {n aS ) = -2\ a lm{(c J b aS )} - 2'y aS (n a s) 
d t (n c ) = - 2 A a (Im{(c 65 )} - Im{(c f 6 a s)}) 

-2q c ((n c ) - n 0 ) 

d t {cb s ) = -iX a ((n s ) + (n c ) + (b s b a s )) 

-(7c + 7s + iu 0 ){cb s ) - iX a 
d t (c ] b aS ) = -iX a ((b s b aS ) ~ (n c ) + (n a s)) 

-('Yc + laS+ iu 0 )(c t b aS ) 

d t (bsKs) = -iK({^b a s) + ( cb s )) 

-( 7 s + 7aS + i2u 0 ) (b s b a s). (11) 


where A Q = A|a|. Initially we will be setting A 5 = A a s = A and 7 s = 7 aS = 7 - The dynamics 
will depend on how the incident held is implemented. In case of a pulsed laser, it is enough 
to redefine the coupling strength as a time-dependent function |Aq,( t )| 2 = X 2 \a\ 2 f(t). f(t ) 
can, for instance, be modelled as a temporal gaussian profile, such as /(f) = exp(— t 2 /2a 2 ), 
where the coherence time is proportional to cr. In that case there isn’t a steady state and, in 
order to analyze the behavior of the held population as a function of the laser power P^, we 
might consider measures at the excitation time, that is, the time at which the population 
gets its hrst maximum, or the time-averaged values of the observables dehned in Eq. as 
(0) t = liniAt-voo So\0) t dt/lXt. 

We hrst focus on the case of a continuum pumping, i.e., A a (t) = X a , on the material. In 
this case, the system of differential equations 0 has a steady state solution due to the 


environment-induced relaxation. To find them, we write the system in (11) as a vectorial 


diffential equation d t x = Mx + b, where M is a time-independent square matrix and 

T T = ((n 5 ), (fias), (h c ), (cbs), (c r b j s ), (c ] b aS ), (cb[ s ), (b s b aS ), (blb[ s )) 
b r = (0, 0, 27 c n 0 , —iX a ,iX a , 0, 0, 0, 0), 


( 12 ) 

(13) 
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and set dtx ss = 0, where ss stands for steady state. To solve this, we invert the matrix M, 
which is only possible if M is not singular, i.e. if its determinant Det(M ) is different from 
zero. The solutions are given by x ss = —M~ 1 b. By simple inspection, it is straightforwardly 
noticed that the Stokes, anti-Stokes and phonon populations satisfy an equilibrium condition 


7 

(hc)ss Uo (y^S/ss (fias')ss) ■ 

7c 


(14) 


This equation shows that if Stokes and anti-Stokes are produced at the same rate, the only 
phonons left are the thermal ones. On the other hand, since in general Stokes production 
is favoured over anti-Stokes (as it becomes clear in a few lines), whenever these scattering 
processes take place, the distribution of phonons in the sample is not determined solely by 
temperature anymore. 

One of the central results in the paper is to show the behaviour of the Stokes-anti-Stokes 
held intensities in terms of the power of the pumping laser. In the steady state regime, we 
can obtain an analytical expression for those intensities, Is( a s), which are proportional to 
the average population of the photon modes (hs(aS))ss- For S- and aS-modes we find 

2(1 + n 0 ) 7 (l + 7 )A 2 — (3 + 4n 0 + 2 y)A „ 


(ftS ) ss 


2(7 + 7 2 - 2A2)(y(l + 7 ) - (1 + 27 )^) ’ 


(15) 


2 wq 7(1 + 7) A 2 - (1 + 4 n 0 - 27) A* 

2(7 + 7 2 - 2Aq)(7(1 + 7) - (1 + 27)^) ’ 


(16) 


where 7 = 7 / 7 c and A Q = Note that A^ = ( X 2 /A)Pl■ The above solutions are well 

behaved, i.e., Det(M) 7 ^ 0, and are physically meaningful if (ns( a s))ss > 0; these conditions 
are fulfilled if \ a < ^ 7(7 + 1)/(1 + 27 ). Therefore our description is valid below the laser 
power upper bound 


Po = A 


7(7 + 1) 
A 2 (l + 27 ) ’ 


(17) 


Note that for low laser power both {ns( a s))ss increase linearly with Pj. In order to easily 
identify different regimes for Pl, the intensities ratio I a s/Is is a good figure of merit. If 7 S> 
{no, 1 }, a condition fulfilled in the analysed experiments, and X a ^ 0 , the ratio (n as ) ss /(n s ) ss 
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reduces to 


{n s ) 


n o 


n 0 + 1 

n 0 
n 0 + 1 

np 

n o + 1 


X 


1 + -2-fk 

2ng Pq 


1 


Pi 


X 


1 + 


2(no+l) Pq 

1 Pl 


2'Mo -Po 


1 - 


Pr 


2(rio + 1) Pq 


-l 


1 + - + 




1 Pi 


n o {no + 1)/ 2 P 0 no(no + 1) \2 P ( 



where Pq ~ ^.77 C /2A 2 . Therefore, the intensities ratio takes the form 


IaS 


n 0 


- L ao _ 

Is no + l 


1 + —+ 


no {no + 1) 


CsclsPl + 


n‘i 
°SaS 


n 0 (n 0 + 1) 


P't 


(18) 

(19) 

( 20 ) 


( 21 ) 



Figure 3: The figure shows the experimental measurements for the three samples investigated in 
the paper (symbols). The lines are the equation (7) in paper plotted using the fitting parameter 


of the table. Dashed black line corresponds to Eq. (7), red line to Eq. (21) only considering the 


quadratic part. Light-blue line correspond to the full equation Eq. (21). 


with CsaS = l/2Po = A 2 /v 4.77 c . We chose to not use this notation in the paper to keep the 
equations related to the fundamental constants. Here C = C a s/Cs, with Cs( a S) defined as 
the proportionality constant such that Is(aS) — Cs(aS){ns(aS))ss- 
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In figure [3] we show that the quadratic term in Eq. (21) does not introduce strong cor¬ 
rections. In fact, all the experiments are done in the limit Pl/Pq 1 so it is reasonable 
that the linear term in this parameter dominates. The blue-dashed lines corresponds to the 


full equation ( 21 ) and the black-dashed line only to the linear part of it, both using the 
fitting parameters of the table presented at the end of the paper. It is clearly seen that these 
both lines basically coincide, further justifying the approximation taken in Eq. (7) of the 
manuscript. 

The phonon population can also be analytically obtained from the stationary analysis 
and it is given by 


= n 0 7(l + 7) + (7 - n 0 )(X 2 /A)P L _ ^ 1 + 27777 
7(1+7)-(1 + 2t)(A VA)P l ~"° 1-a 


( 22 ) 


which is only valid if X a < \/ 7(7 + 1)/(1 + 2y). In the case 7 1, the phonon population 

reduces to 


(n c ) ~ no( 1 + ( I- 2 )Cso.sPl H C 2 s as pj\ (23) 

V n o n o J 

Note that if Pl 7 ^ 0, which activates the S-aS modes, the phonon population is always larger 
than no, therefore the difference ( hs) S s — ( h a s)ss > 0 . 

There are additional effects that allow us to sophisticate our effective description. For 
example, we can also consider effects of phase noise caused by elastic scattering of generated 
phonons with different momenta. This is done by adding an extra Lindbladian term 


C'c(p) = - 2 h c ph c ). 


(24) 


It is easily proven, using Eq. (10), that phase noise does not modify the equations for the 


populations in Eqs. 15 and 16, therefore, the equilibrium condition is preserved. Its main 
changes are for equations of the expected values of the type ( cbs ) and (c^b a s), introducing 
expectation values of second-order operators, as for instance ( cbsh c ). The new system of 
equations is not trivial and analytical solutions are no longer easy to obtain. In that case, 
and for the computation of correlation functions (byr x ,b x "b x iii) (x = S,aS), we compute 
the stationary solution for the density operator p ss , for which dtp ss = 0. We must thus 
solve AotaiPss = 0, which requires an expansion in the Fock basis {\ns, n a s, n c )} G H = 
1~L S ® T-L aS ® T-L c . This basis increases exponentially and must be truncated for practical 
implementations. This is justihable by the fact that since Pl/Pq is always much smaller 
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than one, the steady state photonic population is never large. After numerically obtaining 
p ss we can compute expected values as (0) ss = T r(Op ss ), and we are able to analyze 
changes on the population of the photonic modes and S-aS correlation functions in terms of 
the system parameters. 

MORE THEORETICAL RESULTS 

Figure [4] presents how populations depend on the laser power Pl for 7 = 100 and different 
thermal phonon numbers no. For n$ 7^ 0 and low laser power, both photonic modes increase 
linearly with Pi , with a slope proportional to (no + 1) for S-modes and no for aS-modes. 
In the limit of extremely low temperature (n0 —> 0) the aS-modes always increase with Pf. 

But for a non-zero temperature there are two regimes for the growth of the aS population, 
i.e., I a s oc P L for low laser power and I a s oc Pf for high laser power. Note that for Pl ^ 0 
there exists activation of S- and aS'-modes due to the inelastic scattering, however random 
generation of photonic modes is also expected due to the presence of the thermal bath of 
phonons. The interplay between these two processes will affect the production of photon 
pair S-aS, i.e. the SaS event. 

The two regimes are also reflected in the results for cross-correlation function g 2 ( 0), as 
shown in Fig. [5] Note that the presence of phase noise (7' 7^ 0) does not affect relevantly 
the phenomenology presented in this work (see Fig. [5]). Nevertheless, for very large values 
of Pl the convergence of our numerical fails due to the truncation of the Hilbert space, and 
the results are not conclusive due to the lack of sufficient resources for tackling that regime. 

So far our analysis has been done for continuous pumping, which allows for well behaved 
steady state solutions. When considering a time-dependent but short-time pumping, the 
system response is similar in shape to the implemented pulse and an steady state is not 
expected. In that case, as mentioned before, we computed the long-time averages of the 
observables of interest. The results for the photonic populations are shown in Fig. [6] in 
comparison with the steady state solution for continuum pumping. 

It is clearly seen that the phenomenology studied throughout this paper is preserved, that 
is, the different regimes of the populations with the increasing laser power Pl is preserved, 
and only changes in the intensity of the S- and aS’-modes are observed. The results in Fig. [6] 
allow us to generalize the findings presented along this work. 
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Figure 4: Panel (a) presents Stokes (S'), anti-Stokes ( aS ) and phonon (Ph) populations as a func¬ 
tion of the laser power Pi in units of Po- The different curves correspond to differents temperatures 
set by the thermal average number no (see legend in panel (b)). In panel (b) we show the ratio 
between aS and S mode population as a function of Pi and no- For this figure we set 7 = 100, 
7 ^ = 0 and a continuum pumping. 

We finally want to remark that our theory is based on phenomenological constants that 
will depend on the material under study. The validity of the approximations made to obtain 
analytical expressions for the ratio I a s/Is and the phonon population (n c ) is supported 
by typical and relevant parameters as n 0 , 7 and 7 C . In typical experiments, like the one 
reported in Ref. m, n 0 = 1.6 x 10 3 at T = 295 K, 1 / 7 C « 10 12 s (phonon lifetime) and 
I /7 ~ 130 x 10 -15 s (laser pulse width). The 7 / 7 c is at least one order of magnitude larger 
than one, and four orders of magnitud larger that no- The incident held has a power range 
from 1 mW to 1W and wavelength 785 11111 . The number of incident photons is of the order 
of \a\ 2 ~ 52 x 10 6 — 52 x 10 9 per pulse. 
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Figure 5: (a) Stokes (upper lines) and anti-Stokes (lower lines) modes populations as a function of 
the normalized laser power (Pl/Pq) f° r different values of '/ c /'y c (different line colors, see legend), 
(b) The effect of 'y' c /'y c on g 2 (0) for different values of no (see legends). Other parameters are same 
as those in figure [4] 

ABOUT THE ANALYSIS OF THE EXPERIMENTAL DATA 

A fair question one can ask is wether the experimental results for I a s/Is from diamond [T] 
and graphene [2] cannot be satisfactorily fit considering just the Bose-Einstein distribution 
function. In other words, does the data on the power dependence of I a s/Is from diamond 
and graphene provide convincing case that anything beyond laser heating has been observed? 

The presence of the SaS correlation in diamond is unquestionable, given the results of 
the correlation function g 2 (0 ) [T|. As discussed in our paper, our theory provides perfect 
description for the excitation power dependence of 4s, Is, hs/Is and g 2 (0). The AB-BLG 
is fit in our work considering that the I a s/Is is ruled by thermal effects and, although 
relatively good fitting can be obtained with a small contribution from the SaS phenomenon, 
we raise no controversy in this case. We now discuss results for the tBLG sample [2|, which 
according to our analysis, is dominated by the SaS phenomenon. 

Figure[7]is a plot of what would be the sample temperature T, if extracted directly from 
the Ias/Is experimental data and using the Bose-Einstein phonon distribution function, 
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Figure 6: Stokes (lines) and anti-Stokes (lines with squares) modes populations as a function of 
the normalized laser power ( Pj j jP q ) for different types of incident laser (see legend). The inset 
shows the time profile of the implemented lasers. The parameters are same as those in figure [IJ 
with no = 10~ 3 . 


which from Eq.(l) in the paper is given by: 


T 


ks 


Eph 

In C - In ^ 

is 


(25) 


where E p h is the phonon energy, ks is the Boltzmann constant, and C is the proportionality 
constant that accounts for all the optical properties of the setup and sample. C is chosen 
such that T —» 295 K as Pl —> 0. For tBLG we obtain C tBLG = 26 and for the AB-BLG 
we obtain C AB ~ BLG = 4. The fact that C tBLG > C AB ~ BLG is expected, since tBLG was 
engineered for specific resonance behavior. However, the fact that T is generally larger for 
tBLG, as compared to AB-BLG (see data in Fig. [7]) is not expected. First, the tBLG was 
engineered for resonance with the anti-Stokes photon emission. Therefore it is expected 
that this sample would have a stronger cooling channel, and it would actually exhibit lower 
temperatures than the AB-BLG. Second, this expectation is proven to be correct by the 
lower shift in phonon frequency observed for tBLG, as shown in Ref.[2j. Therefore, the 
tBLG can not be hotter than the AB-BLG, and the result obtained with Eq. ?? is proven 
to be inconsistent. 
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P L (mW) 

Figure 7: Expected effective temperature T for tBLG (red) and AB-BLG (blue), extracted directly 
from the I a s/Is data [2], according to Eqj2l] neglecting the SaS phenomena. 

One can also argue on how sensitive the data fitting is with respect to the parameters 
choice. Parameter values outside the errors given in table I of the paper do not fit the data 
properly. Figure[8] shows three possible fittings of the tBLG data. The red line is the best fit, 
as given in the paper. The black and blue lines were obtained by considering that the Ias/Is 
is ruled by pure SaS phenomena and by pure thermal phenomena, respectively. Notice that 
these two fitting options do not fit the data properly. The worse fitting is obtained using 
the pure thermal phenomena, demonstrating the predominance of the SaS phenomena. 

Figure|9] presents the separate contribution from the thermal phonons (dashed blue lines) 
and from the SaS phenomena (dotted black lines) to the I a s/Is data fitting (solid black lines) 
for diamond (a), tBLG (b) and AB-BLG (c). As discussed in the paper, for diamond and 
tBLG, the SaS phenomena dominates until large laser powers, when the thermal phonons 
overcome. For AB-BLG, the SaS phenomena is minor. 

It remains to be discussed the effect of A5 7^ A a s, which is important when resonance 
effects are in place, i.e. when either Stokes or anti-Stokes has a stronger coupling to the 
pumping field. Considering ^ = e and setting X s = A and A a s = eX, once again as long as 
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Figure 8 : Three different fittings for the I a s/Is data from tBLG [2]. The red points are the 
Ias/Is data. The red line is the same fitting as in the paper. The black line is a fit to the data 
using Eq.(7) from the paper, with C' = 20, (A 2 / 77 c )Al _1 = 1.8 x 10 - 3 mW -1 and T = 295K. The 
blue line is a fit to the data using Eq.(7) from the paper, with C' = 20, (A 2 / 77 c ) 7 l -1 = 0 and 
T = (295 + 35 P L ) K. 


A 7 ^ 0, Eq.(7) in the paper changes to: 

{n a s)ss =£ 2 ( 2 ^ 07(7 + 1 ) - (1 - 2 y + n 0 (3 + e 2 ))A 2 

(' n s)ss y2(n,0 + 1)7(7 + 1) — (1 + n o + (2 + 3 no + 27)s 2 )A 2 

For 7 {no, 1,£ 2 }, this expression is well approximated by 


{^aS ) ss 

(^.s) ss 


£ 2 n 0 
n 0 + 1 





_£1_^ 

7(no+l) J 


xi 


(26) 


which approximately reduces to: 

Mss ^ c 2 n 0 f + 1_ A|\ / + e 2 A 2 \ . 
(ns)ss n 0 + 1 y n 0 7 J y n 0 + 1 7 J ' 

Finally, the intensities ratio, in terms of Pq = A^^ c /2\ 2 and when with e 

given by 


(27) 


A a s/As < 1 , is 


IgS 

is 


C" £ 


/ _2 


no 


- £2 
+ n 0 + (n 0 + 1 ) 


CsaS^L + 


n 0 (n 0 + 1) 


7^2 p2 

^SaS r L 


(28) 


n 0 + 1 
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with CsaS — 1/2-Po- It is straightforward to compute the other case, i.e., when A s < A a s- In 
that case if we redehne e = \ S /X a s < 1, the expression for the intensities ratio is: 



In non-resonant scattering, A $ = A a s is expected due to time-reversal symmetry. This is 
the case for diamond [I]. If resonance with electronic states exists, but there is not much 
difference between the resonance with the S and aS photon emission, then A 5 r\j X a s • This is 
the case of AB-BLG, where the density of states does not vary too much between the S and 
aS photon energies |2j. However, in the case of tBLG, the electronic structure is engineered 
to exhibit specific resonance with the aS photon emission [2], i.e. A a s > A 5. Therefore, 
strictly speaking, Eq. ?? from this Supplemental Information should be used to fit the data, 
rather than Eq. 7 from the paper. The value of £ is easily obtained from the data if we 
consider that this parameter rules the enhancement of C from tBLG as compared to AB- 
BLG. From table I in the paper and Eq. ?? from this Supplemental Information we obtain 
e 2 = 3.5/20. We now add £ = 0.418 in Eq. ?? from this Supplemental Information and 
fit the tBLG data with C' = 3.5, which is consistently the value for AB-BLG. This means 
that the effect of resonance in tBLG is now accounted for with e. The data fitting gives 




10 


10" 2 10" 1 10 ° 10 


10" 2 10" 1 10 ° 10 1 10 2 


P L (mW) 


P L (mW) 


P L (mW) 


Figure 9: Contribution from the thermal phonons and from the SaS phenomena in the fitting 
of Ia,s/Is from (a) diamond, (b) tBLG and (c) AB-BLG. The symbols are the experimental data. 
The solid black lines are the best fits using Eq.(7) from the paper, with the parameters shown in 
Table I. The dashed blue and the dotted black lines give the separate contribution from the thermal 
phonons and from the SaS phenomena, respectively. 
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(A 2 /77 c)^” 1 = (8.1 ± 0.5) x 10 _3 mW _1 and (T[K] = 295 + (15 ± l)P/:[rnW]). Therefore, 
considering the possibility of A s ^ A a s i n our analysis actually enhances the importance of 
the SaS phenomena in tBLG by almost one order of magnitude (compare with the value for 
(A 2 /77 c )7l~ 1 in Table I). Furthermore, the dependence on temperature decreases to half of 
the value obtained for AB-BLG, which is (T[K] = 295 + 30 Pl[ulW]), consistent with the 
G band frequency dependence observed in Fig. 7 of Ref. [2]. Although this picture is fully 
consistent, we chose to keep the simpler equation in the paper since the general picture does 
not change, and the discussion is simpler. 


[1] Kasperczyk, M. Jorio, A. Neu, E. Maletinksy, P. and Novotny, L. Stokes-anti-Stokes Correlations 
in Raman Scattering from Diamond Membranes, http://arxiv.org/abs/1503.00467 

[2] Jorio, A. et al. Optical-Phonon Resonances with Saddle-Point Excitons in Twisted-Bilayer 
Graphene. Nano Lett. 14, 5687 (2014). 


23 




